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This paper presents a combined analysis of the KARMEN and LSND ^ e -carbon cross section 
measurements within the context of a search for u e disappearance at high Am 2 . KARMEN and 
LSND were located at 17.7 m and 29.8 m respectively from the neutrino source, so, the consistency 
of the two measurements, as a function of antineutrino energy, sets strong limits on neutrino oscilla- 
tions. Most of the allowed region from the u e disappearance analysis of the Gallium calibration data 
is excluded at >95% CL and the best fit point is excluded at 3.6cr. Assuming CPT conservation, 
comparisons are also made to the oscillation analyses of reactor antineutrino data. 
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I. INTRODUCTION 

This paper presents an analysis of the i/ e -carbon cross 
section data from the KARMEN [J and LSND [5] experi- 
ments, within the context of electron neutrino oscillations 
at high Am 2 . In a two- neutrino oscillation formalism, the 
probability for v e disappearance is given by: 

P = 1 - sin 2 26» sin 2 (1.27Am 2 (L/£;)), (1) 

where 9 is the mixing angle; Am 2 = to 2 . — m 2 , in eV 2 , 
is the difference between the squared neutrino masses; 
L, in m, is the distance from the neutrino source to the 
detector; and E, in MeV, is the neutrino energy. This 
analysis exploits the fact that KARMEN and LSND were 
located at L =17.7 m and 29.8 m respectively. We use the 
consistency between the cross section measurements to 
place strong constraints on v e disappearance at Am 2 ~ 
1 eV 2 . 

This study is motivated by recent results that can 
be interpreted as oscillations with Am 2 ~ 1 eV 2 . The 
strongest evidence comes from the LSND experiment, 
which observed a — > v e signal corresponding to an 
oscillation probability of (0.264 ± 0.067 ± 0.045)% 0. 
Recent MiniBooNE antineutrino data [4] are in agree- 
ment with LSND, although with less significance, while 
the MiniBooNE neutrino data do not support — > v e 
oscillations [5]. 

High Am 2 muon-to-electron flavor appearance cannot 
be explained in a three neutrino mixing model that also 
incorporates "solar" and "atmospheric" oscillations [BJ. 
As a result, these data have inspired models with three 
active and one sterile (3+1) or 2 sterile (3+2) neutrinos. 
Sterile neutrinos (y s ) do not interact via the weak inter- 
action, but can mix with and cause oscillations between, 
the active flavors. These models predict a v e —> v s signal 
[7] with a large Ato 2 (on the order of a few eV 2 ) com- 
pared to the splittings between the light states (of order 
- 10~ 3 and 10" 4 eV 2 ). Therefore, one can take the three 
light states to be effectively degenerate. 

This degeneracy simplifies the 3+1 model to an ap- 
proximate two-neutrino oscillation model for both ap- 



pearance and disappearance. As a result, Eq. [TJwill be 
applicable to the following discussion, where we will use 
9 ee as the mixing angle relevant to u e disappearance. 

Recently, a reanalysis of reactor v e flux predictions [5] 
has provoked further interest in electron flavor disappear- 
ance in 3+1 models [91411] . This new analysis resulted 
in a shift of the ratio of reactor data-to-prcdiction from 
0.976 ± 0.024 to 0.943 ± 0.023. This deficit with respect 
to prediction is called the "Reactor Anomaly" in this pa- 
per. This can be taken as indication of v e — > D s in a 3+1 
model at 98.6% CL [pj. The best fit is Am 2 = 1.78 eV 2 
and sin 2 26» ee = 0.088 [II]. 

Indications of v e disappearance have arisen from cali- 
bration data taken by the SAGE [12] and GALLEX [13] 
experiments. These used megacurie sources of 51 Cr and 
37 Ar to calibrate the ^ e + 71 Ga — > 71 Ge + e~ experiments. 
The data from SAGE and GALLEX are consistent, and 
show a measured-to-predicted ratio of R = 0.86 ± 0.05 
[TP] . We refer to this as the "Gallium data" in this paper. 
This can be interpreted as a 2.7tr indication of u e — > v s 
oscillations jTOj [14] . The best fit in a 3+1 model corre- 
sponds to a Am 2 = 2.24 eV 2 and sin 2 26 ee = 0.50 [TO] . 
This apparent v e disappearance signal leads to the ar- 
gument |15) that v e —¥ v s , consistent with the Gallium 
data, must be applied to the intrinsic v e background of 
the MiniBooNE — > u e search [S]. 

The results of the 3+1 models have provoked substan- 
tial interest within the community, but they are not de- 
cisive. As a result they are a prime motivation for fur- 
ther studies of electron neutrino disappearance |16[ 117) . 
There are few opportunities for precision v e disappear- 
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FIG. 1: Energy distribution of neutrinos in a DAR beam. 
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ance searches, since most beams have large uncertainties 
in the normalization and energy distribution for the v e 
and v e beam contents. However, decay-at-rest (DAR) 
neutrino beams can provide a unique window on electron- 
flavor neutrino oscillations. 



II. KARMEN AND LSND 

KARMEN and LSND were DAR experiments that ran 
in the 1990's using 800 MeV protons on target. The 
isotropic DAR flux, shown in Fig. [TJ has equal v^, 
and v e content with a well-understood energy spectrum 
described by weak decay physics. The beam energy ex- 
tends to 52.8 MeV. The normalization is known to 10%, 
with the uncertainty dominated by the pion production 
rate per incident proton (THJ, [TI5] . 

KARMEN ran at the ISIS facility at Rutherford Lab- 
oratory, with 200 fiA of protons impinging on a copper, 
tantalum, or uranium target. The center of the nearly 
cubic detector was located at 17.7 m from the proton 
target, at an angle of 100°. The liquid scintillator target 
volume was 56 m 3 and consisted of 512 optically indepen- 
dent modules (17.4 cm x 17.8 cm x 353 cm) wrapped in 
Gadolinium-doped paper. More details are available in 
Ref. [20]. 

LSND used protons from the LAMPF accelerator at 
Los Alamos National Laboratory (LANL), where a 1 mA 
beam of protons impinged on a water target. The center 
of the 8.75 m long, nearly cylindrical detector was lo- 
cated at 29.8 m from the target, at an angle of 12° from 
the proton beam direction. This was an unsegmented 
detector with a fiducial mass of 167 tons of oil (CH2), 
lightly doped with b-PBD scintillator. More details are 
available in Ref. [21]. 

Both experiments measured v e ± 12 C — ^ 12 N gs + e~ 
scattering. In this two-body interaction, with Q-value 
of 17.3 MeV, the neutrino energy can be reconstructed 
by measuring the outgoing visible energy of the electron. 
The 12 N ground state is identified by the subsequent /? 
decay, 12 N gs — s- 12 C + e + + u e , which has a Q-value of 
16.3 MeV and a lifetime of 15.9 ms. 



III. THE KARMEN AND LSND CROSS 
SECTIONS 

The KARMEN and LSND cross section measurements 
for v e + 12 C — > 12 Ng S + e~ [TJI2], m energy bins, are com- 
pared in Fig. [2] The corresponding flux-averaged cross 
sections measured by KARMEN and LSND are given in 
Table |TJ For completeness, we also list the flux-averaged 
cross section for the LANL E225 experiment [22], which 
was located 9 m from a DAR source. (E225 did not 
publish energy-binned cross section measurements.) The 
agreement between all three experiments is excellent. 

The measured cross sections are compared to predic- 
tions by Fukugita, et aL[23] and by Kolbe et al. [2(5] in 
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FIG. 2: The KARMEN (points) and LSND (crosses) mea- 
sured cross sections with statistical errors for v e ± 12 C -V 2 
N ffs + e~ compared to the theoretical prediction of Fukugita, 
et al. (solid line), based on the EPT model, and Kolbe, et al. 
(dashed line), based on the CRPA model. 



Experiment (Dist.) 


Flux-Averaged Cross Section 


Ref. 


KARMEN (17.7 m) 
LSND (29.8 m) 
E225 (9 .0 m) 


(9.1 ±0.5 ±0.8) x 10~ 42 cm 2 
(8.9 ±0.3 ±0.9) x 10~ 42 cm 2 
(1.05 ± 0.10 ± 0.10) x 10~ 41 cm 2 


m 
m 
m 


Prediction 


Flux-Averaged Cross Section 


Ref. 


Fukugita et al. 
Mintz et al. 
Donnelly 
Kolbe et al. 


9.2 x 10~ 42 cm 2 
8.0 x 10 -42 cm 2 
9.4 x 10~ 42 cm 2 
8.9 x 10 -42 cm 2 


m 

M 

21 i 


TABLE I: Top: Flux-averaged u e + 12 C -> e+± 12 N 9S cross sec 



tion measurements with statistical and systematic error. Bot- 
tom: Flux-averaged predictions from EPT (Fukugita, Mintz 
and Donnelly) and CRPA (Kolbe) models. Flux-average cross 
section values are equivalent to those for a neutrino of 35 MeV 
energy. 

Fig. [2] Both models follow an (E v — Q) 2 form, where 
Q = 17.3 MeV. This energy dependence arises because 
the interaction is an allowed transition, converting the + 
( 12 C) state to the 1+ ( 12 N) state. The Fukugita predic- 
tion is calculated within the "elementary particle model" 
(EPT) and has an associated 12% normalization uncer- 
tainty [23,. Other EPT predictions, include Donnelly [23] 
and Mintz, et at, are given in Table [i] For comparison, 
the Kolbe, et al, calculation [26] is performed within a 
"continuum random phase approximation" (CRPA) ap- 
proach. A discussion of the relative merits of EPT ver- 
sus CRPA models for describing this process appears in 
Ref. |27j . From a strictly experimental point of view, 
both EPT and CRPA models fit the data well. To be 
clear, these theoretical results are true predictions rather 
than fits, since they were published well before the KAR- 
MEN and LSND results. 



IV. CONSTRAINTS ON ELECTRON 
NEUTRINO DISAPPEARANCE 

The allowed regions for v e — > v s oscillations are deter- 
mined from a comparison of the LSND and KARMEN 
data with respect to the Fukugita prediction. For a given 
oscillation hypothesis (Am 2 and sin 2 2# ee ), we calculate 
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FIG. 3: Comparisons of the data to various oscillation predic- 
tions for the LSND (top) and KARMEN (bottom) data using 
the Fukugita prediction, as described in the text. 



a combined x 2 for LSND and KARMEN with respect 
to the prediction using the statistical error for each data 
point and employing three pull terms as a method to in- 
corporate systematic uncertainties. The first pull term 
represents the correlated normalization error. As noted 
in the KARMEN paper [J, LSND and KARMEN have 
a 7% systematic error on the neutrino flux normalization 
from the flux simulation that is correlated between the 
two experiments |18[ 1 1 9] . This is combined in quadra- 
ture with the 12% systematic error on the normalization 
for the Fukugita prediction to give the correlated nor- 
malization pull term in the x 2 calculation. The remain- 
ing uncorrelated normalization uncertainties for each ex- 
periment are 7% for LSND [2] and 5% for KARMEN 
PQ. These uncertainties are used as the two other pull 
terms in the x 2 calculation. To determine the 90% CL 
allowed regions in Am 2 and sin 2 2# ee , we marginalize 
over the three normalization pull parameters and use a 
Ax 2 > 4.61 requirement for the two-degrees-of-frecdom 
excluded region. 

The results of the fits using the Fukugita prediction 
are shown in Fig. [3] Table [n] reports the x 2 and degrees 
of freedom (DOF) for various joint fits to the LSND and 
KARMEN data points. The fit without oscillations (No 
Osc) , shown as the long-dashed line in Fig. [3J has a Ax 2 
probability of 91.5% and is only excluded at the 1.7<r 
level. As a result, we use the data to set a limit on v e 
disappearance and calculate the 95% CL exclusion region 
shown in Fig. [4] The best fit, indicated by the solid lines 
in Fig. § is at Am 2 = 7.49 ± 0.39 eV 2 and sin 2 26 ee = 
0.290 ±0.115. 

Comparing the data to an oscillation model with the 
best-fit Gallium parameters illustrates the strong dis- 
agreement, though we note that the Gallium fit had a 



Fit Type 


Am 2 (eV") 


sin 2 26 ee 


X7DOF 


Best Fit 7.49 ± 0.39 
No Osc 

Gallium 2.24 (fixed) 
Reactor Anomaly 1.78 (fixed) 


0.290 ±0.115 
0.000 (fixed) 
0.500 (fixed) 
0.089 (fixed) 


5.5/9 
10.4/11 
34.3/11 
10.2/11 


Kolbe et al. 


Cross Section 






Fit Type 


Am 2 (eV 2 ) 


sin 2 28 ee 


X7D0F 


Best Fit 7.49 ± 0.39 
No Osc 

Gallium 2.24 (fixed) 
Reactor Anomaly 1.78 (fixed) 


0.281 ± 0.115 
0.000 (fixed) 
0.500 (fixed) 
0.089 (fixed) 


6.1/9 
10.7/11 
37.8/11 
10.8/11 



TABLE II: Results of fits using the Fukugita (top) or Kolbe 
(bottom) cross sections for the predicted energy dependence. 
The resultant or assumed Am 2 and sin 2 29 ee values, along 
with the % 2 an< A degrees of freedom (DOF) for the fits are 
shown. 




- 95% CL Limit (Fukugita et al.) 
-95% CL Limit (Kolb et al.) 

- 95% CL Sensitivity 
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FIG. 4: The 95% v e disappearance limit from the Fukugita 
(EPT) fit (solid, black line) compared to the predicted sensi- 
tivity (dotted line). Also shown is the 68% (darker, shaded 
region) and 90% (lighter, shaded region) contours from the 
Gallium experiments. The dashed line is the Kolbe (CRPA) 
fit. 



rather shallow minimum [T0l [14]. The Gallium fit re- 
ported in Table |TT| and shown as the dashed line on Fig. [3] 
is poor. This point has a x 2 probability of less than 
3.2 x 10 -4 and is, therefore, ruled out at 3.6(7. (The 
Ax 2 for this point has a probability of 5.3 x 10~ 7 , which 
corresponds to a 5. Oct exclusion.) Most of the Gallium 
allowed region, indicated at 68% and 90% CL on Fig. [4j 
is excluded at 95% CL by this analysis. 

As discussed above, all models tend to follow a {E v — 
Q) 2 dependence. Nevertheless, small differences between 
the Fukugita (EPT) and Kolbe (CRPA) predictions, 
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-95% CL Limit (LSND and KARMEN) 
-95% CL Limit (KARMEN only) 
- 95% CL Limit (LSND only) 
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sin 2 2B 

FIG. 5: The 95% v e disappearance limit from the combined 
fit (solid, black line) compared to individual fits to KARMEN 
data (dashed) and LSND data (dotted). This is overlaid on 
the Gallium 90% and 68% CL allowed regions. All fits use 
the Fukugita model. 



shown in Fig. [2] allow a test for model dependence. The 
Kolbe fit proceeds in the same way as for the Fukugita 
model. The resulting \ 2 values for the fits are given in 
Table [TTJ. The comparisons of the fits with the data are 
indistinguishable from those shown in Fig. [3] and so are 
not shown here. This leads to the conclusion that there is 
no substantial systematic effect from the energy depen- 
dence of the underlying cross section model. The 95% 
CL exclusion limit from the Kolbe fit is also shown in 
Fig. [4] as the dashed contour and is very similar to the 
Fukugita contour. 

The excluded region in Fig. [3] is better than the ex- 
pected sensitivity region (dotted contour) calculated for 
an underlying null oscillation hypothesis. As a way of 
quantifying this difference, using the Fukugita model fit, 
the Gallium data point is ruled out in a x 2 analysis at 
3.6cr, while the sensitivity would have predicted that, for 
an average experiment with no signal, the Gallium point 
would be ruled out at 2.8a. This strong limit with re- 
spect to the sensitivity is not unlikely; we find that 11% 
of simulated experiments have a high Am 2 95% CL limit 
for sin 2 29 ee at or below the 0.34 95% CL limit of this 
analysis. These studies also show that A^ 2 is a good 
statistic for determining the exclusion regions since 10% 
of the simulated experiments have a A% 2 value for the 
null oscillation hypothesis greater than 4.61 as expected 
for two degrees of freedom. 

The combined fit to KARMEN and LSND is stronger 
than fits to the the individual data sets because the de- 
tectors are at different distances. Requiring the proper 



L as well as E dependence adds an important constraint 
to the oscillation fit. The results of the individual fits are 
compared to the combined fit in Fig. [5] A fit to only the 
KARMEN data yields a best fit x 2 of 2.46 for 4 degrees 
of freedom, with the parameters sin 2 29 = 0.333 ± 0.130 
and Am 2 = 7.54 eV 2 . The x 2 for the null fit was 7.05 for 
6 degrees of freedom. A fit to only the LSND data results 
in a best fit x 2 of 2.27 for 3 degrees of freedom, with the 
parameters sin 2 29 = 0.209 ±0.331 and Am 2 = 3.90 eV 2 . 
The x 2 for the null fit was 3.29 for 5 degrees of freedom. 

Considering fits to these two data sets separately al- 
lows interpretation of certain features in the combined fit. 
We see in Fig. [5] that the KARMEN data dominates at 
high Am 2 because of statistics and that the large vari- 
ations in the limit are driven by this data set. These 
variations appear because the event energy distribution 
spans a limited range. As a result, there are oscillation 
parameters for which the KARMEN distance allows the 
experiment to either be very sensitive or very insensitive 
to disappearance. The LSND data set allows full oscilla- 
tion at Am 2 of 3.9 eV 2 , but KARMEN data does not. At 
low Am 2 , the sensitivity is dominated by the LSND data 
although the combination with KARMEN is significantly 
better. 



V. BROADER INTERPRETATION OF THIS 
CONSTRAINT 

The limit presented here and the Gallium data repre- 
sent the only electron neutrino disappearance results in 
this Ato 2 range. Comparison to other data sets require 
interpretation within models. Specific global analyses are 
beyond the scope of this paper, however, we can consider 
the impact of the limit, in general. 

CPT conservation requires that v e and v e disappear- 
ance should occur at the same rate. Because CPT con- 
servation is embedded in all field theories, comparison of 
electron flavor neutrino and antineutrino disappearance 
is widely regarded as interesting. The relavent antineu- 
trino disappearance data comes from the reactor experi- 
ments. The recently published Reactor Anomaly hints at 
oscillations, and the parameters from the v e disappear- 
ance data [5] can be compared to the z/ e results assuming 
CPT conserving models. 

In Fig. |6j we overlay the Reactor Anomaly allowed 
region at 68% and 90% CL with the 95% CL limit from 
this paper using Fukugita. In this figure, we use a log 
scale for sin 2 29 so that the reactor allowed region is clear. 
One can see that portions of the allowed reactor space 
will be excluded in fits that require CPT conservation. 
However, the Reactor Anomaly best fit parameters give 
a Ax 2 probability of 90.8%, which lies outside of the 95% 
CL excluded region from the cross section analysis limit, 
and one expects portions of the allowed region to survive 
in a global fit. The comparison of the best fit with data, 
shown as the thin line with dots in Fig. |3j is reasonably 
good and similar to the "no oscillation" case, illustrating 
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FIG. 6: The 95% v £ disappearance limit from the Fukugita 
(EPT) fit (solid, black line) and Kolbe (CRPA) fit com- 
pared to the 90% (shaded region) contours from the Reactor 
Anomaly. 

why the KARMEN/LSND data does not speak to the 
entire region of the Reactor Anomaly. This is consistent 
with good the x 2 /DOF for the best fit point of the reactor 
anomaly, reported in Table [ill 

Comparisons to electron-to-muon flavor appearance 
data is much more model dependent than the compar- 
ison of v e to v e disappearance. In this case, one must go 
beyond choosing whether or not to assume CPT conser- 
vation, to include the number of sterile neutrinos (typi- 
cally, one assumes 3+1 or 3+2 models); the mass heirar- 
chy (a 1+3+1 model has recently been published [TT] as 
a variation of the two sterile neutrino case) ; and whether 
or not to include CP violation. 

Inclusion of the appearance data also, necessarily, re- 
quires reference to muon-flavor disappearance as well as 
electron flavor disappearance, since these are all realted. 
As an example, in a 3+1 model: 

sin 2 2(? Me = 4<7 2 4 [/ 2 4 , (2) 
sin 2 2^ = 4<7 2 4 (l-[/ 2 4 ), (3) 
sin 2 20 ee = 4C/ 2 4 (1 - L7 2 4 ); (4) 



where U a i represents the element in the 4-neutrino mix- 
ing matrix. The disappearance limit we present here, 
when interpreted within Eq. |4j substantially limits the 
range of l7 2 4 , and will considerably reduce the space 
of allowed values for matrix elements that can describe 
muon-to-electron flavor appearance in neutrino mode, as 
may be implied by the MiniBooNE low energy excess [5] . 
While the correspondences between experimentally mea- 
sured mixing angles and the underlying matrix elements 
for two-sterile-neutrino models is more complicated, the 
basic point still applies. 



VI. CONCLUSIONS 

This analysis has used the u e + 12 C -> 12 N 9S + e _ cross 
section data from LSND and KARMEN to constrain the 
amount of v e disappearance oscillations at high Am 2 . 
The good agreement between the data sets and with the 
theory, despite different distances of detectors from the 
source, leads to an 95% CL exclusion region which ex- 
tends down to sin 2 29 ee — 0.34 at high Am 2 and to con- 
siderably lower values for some Am 2 ranges. Comparison 
to another underlying cross section model does not sig- 
nificantly change the excluded region. 

The data are in disagreement with the only other u e 
disappearance data set in this Am 2 region, which comes 
from the Gallium experiments. Large portions of the 
allowed region for v e disappearance analysis of the Gal- 
lium calibration data are ruled out at >95% CL. As a 
benchmark, the best fit point for Gallium is excluded at 
3.6cr. This new limit also severely restricts models ad- 
dressing the MiniBooNE results, such as Ref. [15], that 
incorporate v e disappearance to explain the observed en- 
ergy distribution. 

Assuming CPT conservation, this data set can be com- 
pared to the v e disappearance from reactor data. The 
Reactor Anomaly best fit point in this case is within the 
allowed 95% CL contour but some regions of the allowed 
region in Ref. [5] are excluded. 
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